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Abstract. The purpose of this work is to resume investigation of Galactic thin disk dating using nucleocosmochronology with 
Th/Eu stellar abundance ratios, a theme absent from the literature since 1990. [Th/Eu] abundance ratios for a sample of 20 disk 
dwarfs/subgiants of F5 to G8 spectral type with -0.8 < [Fe/H] < +0.3, determined in the first paper of this series, were adopted 
for this analysis. We developed a Galactic chemical evolution model that includes the effect of refuse, which are composed of 
stellar remnants (white dwarfs, neutron stars and black holes) and low-mass stellar formation residues (terrestrial planets, 
comets, etc.), contributing to a better fit to observational constraints. Two Galactic disk ages were estimated, by comparing 
literature data on Th/Eu production and solar abundance ratios to the model (8.7^ * Gyr), and by comparing [Th/Eu] vs. [Fe/H] 
curves from the model to our stellar abundance ratio data ((8.2 ± 1.9) Gyr), yielding the final, average value (8.3 ± 1.8) Gyr. 
This is the first Galactic disk age determined via Th/Eu nucleocosmochronology, and corroborates the most recent white dwarf 
"*^h | ages determined via cooling sequence calculations, which indicate a low age (< 10 Gyr) for the disk. 

Key words. Galaxy: disk - Galaxy: evolution - Stars: late-type - Stars: fundamental parameters - Stars: abundances 

o 

a. 

\ 1 . Introduction tive nuclides. These are compared with the abundances of their 

L* ' , daughter nuclides, or of other nuclides that are created by 

i> • Current estimates of the age of the Galactic thin disk are .., , .■ ^ 

• i-H , & the same or a similar nucleosynthetic process. Depending on 

\,; obtained by dating the oldest open clusters or white dwarfs. , , f , f e ■ . .■ , 

J\ 1 3 & F the half-life of the chosen nuclide, different timescal es can 

C Ages of open clusters are determined by fitting isochrones, bg probed The Th/Eu chronometer> first d by g^J 

and ages of white dwarfs are determined using cooling se- JTTioTJv a t * tU CA . i a- i 

& & & 1 1989), is adequate to assess the age of the Galactic disk, since 

quences. Isochrone and cooling sequence calculations re- 23 2 Th ig a radioactive nuclide with a 14 05 G half . ufe e > 

quire deep knowledge of stellar evolution, are very com- c ., , c .. , £ , • a\ r- ■ 

M r & 3 of the order of magnitude of the age being assessed). Eu is a 

plex, and depend on a large number of physical parameters .. f . , . e ■ , ... , ■ , 

F ' F & satisfactory element tor comparison, because it is produced al- 

known at different levels of uncertainty. Many important as- , , , nnrt , A - . rr : TTf^vvl u *u 

3 most exclusively (97%, according to Burns et al. 2000) by the 
pects of stellar evolution, like the influence of rotation, are , .... . , „ • ■ 

F same nucleosynthetic process that produces all Th, the rapid 

not well known, and may have a strong influence on the neutron . capture proC ess (r-process). 

outcome of isochrone calculations. Furthermore, white dwarf I 1 j 1 

cooling sequences also depend on calculations of degenerate- Since ® n ^ n carried out the first Th ^ u nu " 

matter physics. Nucleocosmochronology is a dating method cleocosmochronological dating of an ultra-metal-poor (UMP) 

that makes use of only a few results of main sequence stellar star ' f literature has b een virtually dominated by these ob- 

, .. a ^ t u e n - a a «■ iects (Truran et al. 2002, and references therein). The reason 
evolution models, therefore allowing a quasi-independent ven- J 1 *■ ** ' 

c t - ft uf a * u • for such interest lies in the strong simplifications that can be 

fication of the afore mentioned techniques. & F 

applied to the dating of UMP stars, which does not require the 



Nucleocosmochronology estimates timescales for astro- 
physical objects and events by using abundances of radioac- 



use of Galactic chemical evolution (GCE) models. Dating of 
disk stars, on the other hand, does require the use of GCE mod- 
Send offprint requests to: E.F. del Peloso els, and is thus a very intricate process. O wing to the com- 

1 All references to the Galactic disk must be regarded, in this work, plexity of the analysis, Ida Silva et all Jl990h is the only work 
as references to the thin disk, unless otherwise specified. ever published for disk stars, presenting preliminary results for 
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only four objects. This work aims at resuming investigation of 
Galactic disk dating using [Th/Eu] abu ndance ratios. 2 

In the first paper of this series Jdel Peloso et alJ [2005. 
Paper I), [Th/Eu] abundance ratios were determined for a sam- 
ple of 20 disk dwarf s/subgiants of F5 to G8 spectral type with 
-0.8 < [Fe/H] < +0.3. In what follows, we employ a GCE 
model developed by us to determine two Galactic disk esti- 
mates: one by comparing Th/Eu production and solar abun- 
dance ratios from the literature with the model, and one by 
comparing the stellar abundance data determined by us with 
[Th/Eu] vs. [Fe/H] curves obtained from the model. 

2. Nucleocosmochronological analysis 

Stars that compose our sample have been formed all along the 
disk lifetime. Th present in them was synthesized in a number 
of stellar generations, the younger stars receiving contributions 
from a greater number of them; consequently, matter synthe- 
sized by each generation has decayed for different amounts of 
time. On that account, a GCE model is indispensable for a cor- 
rect interpretation of the evolution of the [Th/Eu] abundance 
ratios, from which we derive the Galactic disk age. 

2.1. Galactic chemical evolution model 
2.1 .1 . Basic description 

We have develop ed a GCE model based on 

Page l & Tautvaisiend i 19951 PT95) with the inclusion of 
the effect of re fuse, following an a melio rated version of the 
formulation of Rocha-Pint o et alJ ([1994). The inclusion of 
the effects of refuse, which are composed of stellar remnants 
(white dwarfs, neutron stars, and black holes) and low-mass 
stellar formation residues (terrestrial planets, comets, etc. - 
hereafter referred to as simply residues), contribute to a better 
fit of the Fe abundance evolution model to the metallicity 
distribution of the G dwarfs in the solar neighbourhood. 

The mod el of PT95 is composed of m ultiple phases. 
According to iBeers & Sommer-Larsenl Jl995t) . about 30% of 
the stars with [Fe/H] < —1.5 within 1 kpc of the Galactic plane 
belong to the disk, which shows a considerable intersection be- 
tween the disk and halo metallicity distributions. Several other 
works rev eal this intersection (PT 95, and references therein). 
Moreover. IWvse & Gilmordd 19921) argue that the angular mo- 
mentum distribution function of the halo stars makes it highly 
improbable that the gas removed from the halo during its stellar 
formation phase was ever captured by the disk, but more likely 
by the bulge. Hence, disk and halo initiated their formations 
disconnected from each other, both from the primordial, unen- 
riched gas, and thus the disk can be modelled independently 
from the halo. Our model is not rigorously a Galactic chemi- 
cal evolution model, but rather a Galactic disk chemical evolu- 
tion (GDCE) model. Accordingly, we assume a simple model 

2 In this paper we employ the following customary spectroscopic 
notations: absolute abundance logs(A) = log l0 (N A /Nu) + 12.0, and 
abundance ratio [A/B] = log 10 (N A /A r B)star - log i0 (N A /N B )o, where 
N A and N% are the abundances of elements A and B, respectively, in 
atoms cirT 3 . 



(Ivan den Berghl 1 962l Isdimidtl 1 9631) until ~ 0.037g, where Tc 
is the Ga lactic disk age. After this initial phase, we assume the 
model of lClavtonl dl985l) . with a mass infall rate 



where to is a constant that represents the efficiency of interstel- 
lar gas conversion into stars, g is the mass of interstellar gas, 
u — cot is a dimensionless time-like variable, uq is an arbitrary 
parameter, and A: is a small positive integer. We define a time- 
like variable U\, which corresponds to the transition between 
simple and Clayton models. 

We adopted the d elayed production approximation devel- 
oped by PageJ dl989h . which takes into account the delay in 
production of elements that are synthesized mainly in stars of 
slow evolution, like Fe, which is predominantly generated in 
Type la supernovae (with a typical timescale of 1 Gyr). In this 
approximation, it is assumed that the elements whose produc- 
tion is delayed start being ejected at a time A after the start of 
stellar formation. The abundance A,- of element i is comprised 
of two components: a prompt component An and a delayed 
component Aq, so that A, = An + Aa- Some elements, like 
Th, are generated exclusively by prompt processes, and con- 
sequently have a null delayed component (Aq = 0). Others, 
like Fe, are formed by both kinds of processes, and thus have 
non-null prompt and delayed components. 

The original PT95 model was modified as follows. In 
iRocha-Pinto et alJ i 19941) . it is assumed that residues evaporate 
a considerable amount of H and He, retaining metals and dilut- 
ing the interstellar medium (ISM). This dilution effect is equiv- 
alent, mathematically, to a second sour ce of metal-poor infall , 
which is important for GDCE models (Chiapp ini et al.lfl"997h . 
In our calculations, the parameter F represents the dilution of 
the ISM due to residue evaporation, and is defined as 

(Jovian planet and residue formation rate) 

r = — — , (2) 

(stellar formation rate) - (stellar ejection rate) 

where the stellar formation includes the formation of brown 
dwarfs and the ejection accounts for the returned fraction (to 
the ISM) due to stellar death. 

We assumed that residue formation, and hence the effect 
of refuse, becomes efficient when metallicity becomes high 
enough to allow the coagulation of planetesimals ([Fe/H] ~ 
-0.30). The time-like variable u% which corresponds to this 
moment, was set initially at U2/0J = 0.40 Tc, but was later re- 
vised to 0.42 Tc by fitting observational constraints. The final 
structure of our model is made up of five phases: 

1. u <u\ \ simple model, 

2. u\ < u < coA: Clayton model with instantaneous recycling, 

3. to A <u <u \ + coA: Clayton model with delayed production 
contribution from stars born during the simple model phase, 

4. Mi + u)A < u < U2'. Clayton model with delayed production 
contribution from stars born after the start of infall, 

5. u > ii2'. ISM dilution due to effect of refuse. 

A schema depicting these phases is presented in Fig.[Q 
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Fig. 1. Structure of the adopted GDCE model (not to scale). A, 
is the abundance of element i. The indices 1 and 2 represent the 
prompt and delayed contributions, respectively. The indices s, 
c, c', and r are relative to the simple model (s), Clayton model 
with delayed production contribution from stars born during 
the simple model phase (c) and after the start of infall (c'), and 
ISM dilution due to effect of refuse (r). 

Eu evolution is modeled by the same formulation used for 
Fe. Th requires specific formulation due to its radioactive de- 
cay. Another important phenomenon, besides radioactive de- 
cay, must be taken into consideration when modeling Th evo- 
lution: selectiv e destruction by pho tonuclear reactions in the 
stellar interior ( Malanev et al. 1989). During H burning via the 
CN cycle, 7.29 MeV and 7.55 MeV y photons are emitted. 
These photons have energies below the photonuclear thresh- 
old of most elements existing in stellar interiors, including Fe 
and Eu. However, Th has a threshold of ~ 6 MeV, which al- 
lows the destruction of up to ~ 40% of its initial content in a 
stellar core. When the star evolves to become a red giant, the 
first dredge-up takes Th-depleted matter from the stellar inte- 
rior to the photosphere. Matter ejected at the end of stellar evo- 
lution has a reduction in Th abundanc e of ~ 10%-20%. This 
pheno menon is properly c onsidered inlAranv-Prado & Maciell 
Jl998l APM981. following lMalanev et alJ dl989h . The general 
equation for the abundance A,- is 



d(gA,) 

du 
where 

Ji = 



+ + J 2 )gAi = 



Pi 



1, 

i + r, 
o. 

i+A, 
P\g(u), 



for u < U2, 
for u > «2> 

for stable elements, 

for r-process radionuclides, 



. P\g{u) + p' 2 g(u - wA), 



for u < wA, or for 
prompt-enrichment 
elements, and 
for u > a>A. 



(3) 



(4) 



(5) 



(6) 



Ai is the decay constant and A, is the destruction term (APM98, 
Eq. 9) of element i; p\ and p' 2 are the net yields of element i, 
as defined by Eq. 8 of APM98 (as well as by Eqs. 6 and 7 
of PT95); the indices 1 and 2 refer to the prompt and delayed 
contributions, respectively. 

2.1.2. Observational constraints and adopted 
constants 

Calculations were carried out for four Galactic disk ages: 6, 9, 
12, and 15 Gyr. Parameter u\ was kept at the value assumed 
by PT95 (0.14). The delay A was fixed at 1.1 Gyr both for Fe 



and for Eu; this value is consistent with the average evolution 
timescales of progenitors of Type la supernovae (which are the 
main astrophysical sites of Fe production) and asymptotic giant 
branch stars (which are the sites of s-process Eu production). 
Note that delayed production of Eu is properly accounted for 
in this work, but that its effects are nonetheless s mall, since Eu 
is majorly (97%, according to iBurris et"al]l2000l) synthesized 
by the r-process. 

We adopted the initial mass function (IMF) of iRroupal 



( 2001, Eq. 6), which covers non-stellar masses (brown dwarfs) 
down to 0.010 m Q . We extended this lower segment down to 
0.001 m Q , keeping the same slope, and included a new segment 
corresponding to the residues that evaporate, with a slope j 
which is a free model parameter. The complete IMF was de- 
fined as follows: 



O(m) 



' a\ m 


-j 




-0.3 


02 m 


A3 m 


-1.8 




-2.7 


«4 m 




-2.3 


K «5 tri 



for 0.000 < m/m Q < 0.001; 
for 0.001 < m/m G < 0.080; 
for 0.080 < m/m G < 0.500; 
for 0.500 < m/m e < 1.000; 
for 1.000 < m/m < 62.000. 



(7) 



It is worth noting that the IMF for the residues relates to 
the mass of residues prior to evaporation. If the IMF from 
iMiller & Scald Jl9 79) had been used, the slope of the residues 
mass range would have to be steeper to compensate for the de- 
crease at low mass range, maintaining normalisation. The sud- 
den change of slope between residue and brown dwarf ranges 
would be harder to explain than the slowe r, more continuous 
increase that can be achieved when using Kroupa's IMF. The 
same would occur if we adopted the IMF from Salpeter ( 1955). 

The parameters u> and uq were evaluated for different val- 
ues of k and for each Galactic disk age, so that at Tq we 
have p = present-day gas fraction = 0.11 and g = mass of gas 
~ (0.73 initial disk mass), as considered by PT95. For a given 
set of values of the parameters k, u, u , and Tq, we fitted the 
metallicity distribution of G dw arfs obtained from the model to 
the observational distribution o flRocha-Pinto & Maciell 11996). 
The fitting was accomplished by changing the parameter j 
and, consequently, the parameter T, since the Jovian planet 
and residue formation rate is given by J m<t>(m) dm (cf. 
Equations |3 and 0. For a given value of parameter j, param- 
eters a\, fl2> fl 3> a 4, an d as (Equation^) were obtained by en- 
forcing IMF continuity and normalisation. The best values of 
j (2.562-2.605) are similar to the IMF slopes for the ranges of 
stellar masses that would form planets (2.7 for 0.5 < m/m Q < 
1.0 and 2.3 for 1.0 < m/m G < 62.0), which implies that equal 
masses form stars and residues - see Table |2] In this table, 
f\ = Pi/A Q and f 2 = p 2 /A @ . 

Good fits were found for k — 2 and k = 3; we adopted 
k — 2, because this value provides a better fit to the four 
Galactic disk ages simultaneously. This lends support to the 
hypothesis that the effect of refuse can substitute for the dilu- 
tion caused by infall at later stages of Galactic evolution, since 
PT95 manage to obtain good fits only with k — 3. The final 
metallicity distributions of G dwarfs obtained from the model 
are compared to Rocha-Pinto & M aciell dl996l) in Fig.|2] 

For the modeling of Eu and Th evolutions, the relevant 
parameters are their respective f\ and f%. For Eu, the to- 
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Observational distribution: 
Rocha-Pinto & Maciel 1996 
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r = 0.301 ;T G = 6Gyr 
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r = 0.199; T a =15Gyr 




Fig. 2. Metallicity distribution of G dwarfs in the solar 
neighbourhood, where AN/N is the normalised number of 
stars. The thick curve is th e observational distribution of 
iRocha-Pinto & Macielldl996l) . The thin curves were calculated 
using our GDCE model for the four different Galactic disk ages 
shown in the legend. 



tal value / to taiEu, which takes prompt and delayed production 
into consideration, was determined by fitting the [Eu/Fe] vs. 
[Fe/H] curve obtained from the calculations to our stellar data 
(Table 14 of Paper I). We assumed that prompt production of 
Eu is accomplished by the r-process, while delayed production 
is carried out by the s-process (with relative contributions 0.97 
and 0.03, respectively, according to B urris et alJ2000l) . leading 
to /ieu = 0.97 / tota i Eu and / 2Eu = 0.03 /total e u - Th is synthe- 
sized exclusively by the r-process and its production was con- 
sidered proportional to the production of Eu: /xh = X /total Eu 
and /2Th = 0. The parameter^ used was the one that provided 
the best simultaneous fit of the [Th/H], [Th/Fe], and [Th/Eu] 
vs. [Fe/H] curves, obtained by the model, to our stellar data. 
Table provides summarised descriptions of the main GDCE 
model parameters, and the final values adopted are listed in 
Table|5J as a function of Galactic disk age. 



2.2. The age of the Galactic disk 

We have estimated the age of the Galactic disk in two different 
ways: using our GDCE model along with literature production 
ratio data, and with our own stellar abundance data. 



2.2.1 . T G from literature data 

Inspection of Table [2] shows that parameter x is n °t constant, 
but dependent on the Galactic disk age. Therefore, the age can 
be estimated by constructing &x vs. Tq diagram, fitting a curve 
to the diagram to determine a relation between the two param- 
eters, and then using literature determinations of x- Parameter 
X can be calculated from literature determinations of solar Eu 
and Th abundances and Th/Eu production ratio through the re- 
lation x = 0? Th /.P Eu )/(Th/Eu) Q , where (Th/Eu) is the solar 
abundance ratio today. 



Table 1. Summarised descriptions of the main GDCE model 
parameters. 

Parameter Short description 

k Small positive integer used to parametrize the mass 

infall rate F. 

A Time after the start of stellar formation, after which 

delayed-production elements start being ejected. 
u) Efficiency of interstellar gas conversion into stars. 

«o Arbitrary parameter in the mass infall rate F. 

u\ Time of transition between simple and Clayton 

model, multiplied by w. 
U2 Time when residue formation becomes efficient, 

multiplied by lj. 
u Tg Age of the Galactic disk, multiplied by u>. 
j IMF slope corresponding to residues that evaporate. 

T Dilution of ISM due to residue evaporation. 

fx Fc Fe yield divided by Fe solar abundance today 

(prompt contribution). 
/2Fe Fe yield divided by Fe solar abundance today 
(delayed contribution), 
/total Eu Eu yield divided by Eu solar abundance today. 
fi eu Prompt contribution of /, ot ai En ■ 
f 2 eu Delayed contribution of / totaI E „ . 
fi Th Th yield divided by Th solar abundance today. 
X Relation between / 1Th and / lola i E „. 

Table 2. Adopted GDCE model parameters. For a short de- 
scription of the entries, see Tabled details can be found in the 
text. 



To 



Parameter 6 Gyr 9 Gyr 12 Gyr 15 Gyr 



k 

A (Gyr) 
u> 

I/O 

"i 

i 
r 

/lFe 
J2Fe 
/otalEu 
/lEu 
7*2 Eu 
/iTh 

X 



0.746 
2.040 



1.880 
4.476 
2.605 
0.301 



0.630 
0.611 
0.019 
0.740 
1.175 



2 

1.1 

0.488 0.362 
1.500 1.250 

0.140 

1.845 1.824 
4.392 4.344 
2.586 2.571 
0.251 0.217 

0.240 

0.467 

0.670 0.710 



0.285 
1.110 



1.796 
4.275 
2.562 
0.199 



0.650 
0.020 
0.905 
1.350 



0.689 
0.021 
1.186 
1.670 



0.750 
0.728 
0.022 
1.425 
1.900 



Observations: Ut d = 

/2Eu 



W 7"g; /lEu - 0.97 /total Eu! 
= 0.03 / tot a| Eu J /l Th = X /total Eu • 



From th e met eoritic determinations originally presented 
bv lLoddersI J2003I) we have loge(Th) = +0.06 ± 0.04 and 
logg(Eu) = +0.49+0.04, which leads to (Th/Eu) Q = +0.372+ 
0.069. These abundance s were not taken direct ly from lLoddersI 
(120031) . but rather from Asnlun det al.l (12005). who corrected 
their zero-point in order to make them compatible with their 
Si photospheric value. It is much more difficult to estimate the 
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production ratio p / p . Some of the works from the literature 
that estimate the ages of UMP stars through Th/Eu nucleocos- 
mochronology use the solar Th/Eu abundance ratio, corrected 
for the Th decay since the formation of the Sun, as an estimate 
of the production ratio. This approach would lead to a good 
estimate only if all solar Th were produced in a single (or at 
most a few) nucleosynthesis episodes just prior to the forma- 
tion of the Solar System. This hypothesis is known to be false, 
and stellar ages estimated with production ratios obtained this 
way yield lower limits only. For this reason, it is preferable to 
use theoretical estimates. 

Theoretical production ratios are derived from models of 
r-process nucleosynthesis, and subject to many uncertainties, 
mainly because the astrophysical site of production by this 
process is not well known yet, and because it is not possible 
to obtain laboratory data on many very-neutron-rich elements 
which are very far from the /^-stability valley. A procedure fre- 
quently employed to reduce these uncertainties is to constrain 
the model results to solar abundances. Heavy elements cur- 
rently present in the Sun have been synthesized by both the 
r-process and s-process (with the exception of U and Th, pro- 
duced exclusively by the r-process). R-process nucleosynthesis 
models are then forced to reproduce just the fraction of the so- 
lar abundances that was produced by the r-process. This con- 
straint is valid only if all stars produce r-process elements in 
solar proportions. The abundances of heavy elements in UMP 
stars could also be used to constrain the r-process nucleosyn- 
thesis models, but their high uncertainty makes it preferable to 
use the solar abundances. 

Ma ny work s from the li te rature JSneden et alJ 1 19961 

19981 kOOOalbl IWestinetal] 1200(1 iBurris et alJ [2000; 
Johnson & Boltel 1200 ll ICowan et alJ 1 19991 120021 among 
others) corroborate the so-called universality of the r-process 
abundances. According to this hypothesis, r-process nucle- 
osynthesis would occur at only one astrophysical site and 
would always produce elements in solar proportions. These 
works support the universality hypothesis by the successful 
comparison of heavy element abundances in UMP stars (rep- 
resenting the product of nucleosynthesis from few preceding 
stars) with r-process fraction of solar abundances (representing 
the cumulative production during several Gyr of Galactic 
evolution). This comparison is carried out for the second and 
third r-process peaks only (56 < Z < 72 and 73 < Z < 82, 
respectively). 3 To enable the use of r-process nucleosynthesis 
models which are constrained to solar abundances in the 
calculation of Th/Eu production ratios, Th must also be 
included in the universality. It is usually assumed that, if the 
universality is valid for second and third peaks, it can be 
extended to the actinides. Determination of more precise Pb 
and Bi abundances in the Sun and UMP stars would allow us 



3 Some of the works cited above show that there is no agreement 
for elements of the first r-process peak (Z < 55). This hints at the exis- 
tence of two different sites of r-process nucleosynthesis: one that pro- 
duces low mass elements (first peak), and one that produces high mass 
elements (second and third peaks, and actinides). The non-universality 
of low mass elements does not represent a problem for Th/Eu nucleo- 
cosmochronology, since both Eu and Th have atomic numbers Z > 55. 



to better investigate if this extension is indeed valid or not. 
Since these elements are produced mainly (more than 80%) by 
the ff-decay of elements with atomic mass 209 < A < 255, it 
would be very unlikely that solar Pb and Bi abundances would 
be compatible with those from UMP stars if universality was 
not valid for elements with A > 209. 

The supposition that universality can be extended to the ac- 
tinides is still a matter of strong polemic in the literature. In 
particular, th e analyses of the ha lo U MP giant CS 3 1 082-001 
carried out bv lCavrel et all ll200lh and lffill et all J2002h provide 
data against it. This star presents a solar abundance pattern for 
elements of the second peak, but its third peak does not match 
the Sun. While Os and Ir seem to be overabundant relative to 
the solar r-process fraction, Pb appears to be deficient. A more 
significant discrepancy is that the star has an abundance ratio 
Th/Eu = +0.603, which is much larger than the solar value 
corrected for Th decay (Th/Eu) Q = +0.434; UMP stars are ex- 
pected to present Th/Eu abundance ratios considerably lower 
than the solar corrected value, due to their higher age. This high 
value indicates that the star could have been formed with an in- 
herent actinide superabundance relative to the elements of the 
second peak, which means that universality would not hold for 
Th. CS 3 1 082-001 is the only star hitherto studied that presents 
important discrepancies for elements with Z > 55, whereas evi- 
dence favorable for the validity of universality for these atomic 
masses is becoming increasingly numerous. It is not clear if 
CS 31 082-001 is merely chemically peculiar or if its discrep- 
ancies could be present in other yet unobserved stars; if the lat- 
ter is true, it may not be possible to extend universality to the 
actinides. If Eu and Th are produced in different astrophysical 
sites, the Th/Eu production ratio may have suffered variations 
during the evolution of the Galaxy, provoked by changes in the 
relative frequencies of their respective progenitors. 

In this paper, we assume that universality can be extended 
to Th, considering that indications to the contrary are very 
sparse. But even if constraining the nucleosynthesis models to 
solar abundances reduces the uncertainties, they nevertheless 
remain very high, because the calculations can be based on dif- 
ferent far-from-stability mass models and other th eoretical sup- 
posit ions (e.g., neutron densities'). ICowan et alJ dl997l fl999) 
and ISchatz et alJ J2£)02) derive values between 0.47 and 0.55. 
We adopted the center of the interval, p Th /p Eu = 0.51 ± 0.04, 
arriving at x — 1 .37 1 lS*3«- In Fig. |3 the Tc vs. x diagram 
constructed for our GDCE model is presented. A linear fit to 
the data was used to estimate a Galactic disk age of Tq = 
8.7*3? ^y 1 ' Taking into account the linear fit scatter (cr = 0.5), 
we arrive at a final value Tq = 8.7*|? Gyr. The high uncer- 
tainty of this result reflects the difficulties of estimating the 
Th/Eu production and solar abundance ratios, as well as the 
low sensitivity of the Galactic disk age to parameter ;f. 

2.2.2. T G from our stellar data 

[Th/Eu] vs. [Fe/H] curves were constructed from our GDCE 
model for each Galactic disk age. These curves, along with our 
stellar abundance ratio data, can be seen in Fig. |4] Note that 
the curves draw further apart at high metallicities. This has a 
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Fig. 3. Tc vs. x diagram for our GDCE model. The dotted curve 
is a linear fit. The solid lines display the value adopted from the 
literature for the x parameter (1.371^^) and its respective 
disk age (8.7+5-3 Gyr). 



very important consequence for the analysis: contrary to the 
intuitive notion that the most metal-poor stars are the most im- 
portant, it is the most metal-rich ones which better discriminate 
between different ages. In order to determine the age that best 
fits our data, we calculated a quadratic deviation between our 
data points and each curve, according to the equation 



Total deviation = ^ {[Th/Eu]; - /([Fe/H],)} 2 



(8) 



where the index i represents the z'-th star, and /([Fe/H],) is a 
4 th order polynomial fit to the GDCE curve. We were not able 
to use all 21 sample stars to calculate the deviations because 
the two most metal-rich objects (HD 128 620 and HD 160691) 
fall out of the interval where the curves are defined; that is why 
the total deviations were calculated by summing 19 individual 
stellar deviations, and not 21. 

Once a deviation was determined for each curve, we traced 
a deviation vs. Tq diagram, and fitted a 2 nd order polynomial 
to it (Fig. [5}- The Galactic disk age that best fits our stellar 
[Th/Eu] abundance ratio data was obtained by minimising the 
fitted polynomial: Tq — 8.2 Gyr. 

We estimated the Galactic disk age uncertainty related to 
the abundance ratio uncertainties by Monte Carlo simulation. 
We developed a code that adds Gaussian random errors to the 
stellar [Th/Eu] and [Fe/H] abundance ratios, where the width of 
the Gaussians are the uncertainties adopted for the respective 
abundance ratio. After this, the code calculates the new total 
deviations between the GDCE curves and the modified stellar 
data, fits a 2 nd order polynomial to the [deviation, Tq] data, and 
calculates the polynomial minimum, obtaining a new age. It is 
possible to choose the number of simulations, each time with 
new abundance ratios, added to different Gaussian errors. The 
ages obtained are saved in a file for further analysis. 

We carried out 75 million simulations. The Galactic disk 
age distribution obtained is presented in Fig. |S] where ages 
were counted in 0.5 Gyr bins. The few negative ages obtained 




-0.4 -0.2 

[Fe/H] 

Fig. 4. [Th/Eu] vs. [Fe/H] diagram for all sample stars. Curves 
were calculated using our GDCE model for the four different 
Galactic disk ages shown in the legend. 
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Fig. 5. Determination of Galactic disk age that best fits our stel- 
lar [Th/Eu] abundance ratio data. Points are the total deviations 
between each curve and the stellar data presented in Fig.|4] cal- 
culated according to Equation[8] The curve is a 2 nd order poly- 
nomial fit. The vertical line marks the polynomial minimum, 
corresponding to best Galactic disk age Tq = 8.2 Gyr. 



(< 0.2%) were removed, since they make no physical sense. 
The distribution has a clear increase in noise as the age in- 
creases (after the mode). This happens because the distribution 
was constructed by counting the number of ages between two 
values. Thus, as the age increases, the number of counts de- 
creases, and statistical fluctuation increases. Because of this, 
we truncated the distribution at 125 Gyr, where the counts drop 
below 100 (corresponding to a 10% statistical uncertainty). 
One other issue, besides the noise, forced us to truncate the 
distribution at a point where the counts are still statistically sig- 
nificant. A continuous probability distribution must fall asymp- 
totically to zero. But this cannot occur with our distribution, as 
it was constructed by counting, and counts have no fractional 
values between and 1. Hence, the result of the calculations 
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T G (Gyr) 

Fig. 6. Galactic disk age distribution obtained through a Monte 
Carlo simulation of the [Th/Eu] and [Fe/H] abundance ratio 
uncertainties. Simulation results have been counted in 0.5 Gyr 
bins, and N is the total count per bin. The vertical line marks the 
age at which the distribution was cut (125 Gyr), corresponding 
to a count fall below 100 (i.e., 10% uncertainty). 

performed to derive an age uncertainty would be incorrect if 
we took into consideration the complete distribution. 

A good estimate for our age uncertainty is the width of the 
distribution. There is more than one way to estimate this width. 
The most widely used is the variance, determined through the 
equation 

o 2 ^ —— (r G -r G ) 2 dr G , (9) 

JO JVtotal 

where iV to tai is the total number of simulations within the inte- 
gration limits and r G is the average age. The integration limits 
should have been —cm and +00, but since the distribution was 
truncated at and 125, we carried out the integration accord- 
ingly. In Fig.Qthe integrand of Equation|5]is presented, and it 
can be seen that its right wing falls very slowly. Even if we ex- 
tended the upper integration limit to +00, it is possible that con- 
vergence could not be achieved. Even if we managed to achieve 
convergence, the value obtained would be very large, and not 
representative of the width of the distribution. 

A more robust esti mator of the wid th is the mean absolute 
deviation (MAbD, see lPress et alJl9 92): 

r 125 n 

MAbD = — - |r G -r G |dr G . (10) 

J() total 

As can be seen in Fig.0 the integrand of Equation [TO] falls to 
zero much faster than that of Equation |5] Therefore, the inte- 
gration converges, and the error introduced by truncating the 
distribution can be neglected. Carrying out a numeric integra- 
tion, we arrive at MAbD = 1.9 Gyr. The procedure was re- 
peated several times, with different bin sizes, and we verified 
that the MAbD is not dependent on the bin size. 

The final value obtained using our stellar data is T G = (8.2+ 
1 .9) Gyr. It must be noted that the cited uncertainty is related 
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Fig. 7. Integrands of Equations l9landlT0l used to calculate the 
variance and the mean absolute deviation (MAbD) of the age 
distribution, respectively. 

to the stellar [Th/Eu] and [Fe/H] abundance ratio uncertainties 
alone, and does not take into consideration the uncertainties of 
the GDCE model itself, which are very difficult to estimate. 
The uncertainty related to the model could very well be the 
main source of age uncertainty. 

2.2.3. Adopted Galactic thin disk age 

The adopted Galactic disk age was calculated by combining 
the two estimates obtained with literature data and with our 
stellar abundance ratios. These values were combined using 
the maximum likelihood method, assuming that each one fol- 
lows a Gaussian probability distribution, which results in a 
weighted average using the reciprocal of the square uncertain- 
ties as weights. The final, adopted Galactic disk age is 

FINAL T G = (8.3 ± 1.8) Gyr. 
3. Conclusions 

[Th/Eu] vs. [Fe/H] curves have been constructed for four 
Galactic disk ages (6, 9, 12, and 15 Gyr) from a GDCE model 
developed by us. These curves were compared to the stellar 
abundance ratio data obtained in Paper I in order to determine 
the Galactic disk age. The age that best fits our data was ob- 
tained by minimising the total quadratic deviation between the 
stellar data and the theoretical curves. A Monte Carlo simula- 
tion was carried out to estimate the age uncertainty related to 
the abundance ratio uncertainties. The value obtained was Tq = 
(8.2+1.9) Gyr. The age was also estimated using literature data 
along with our GDCE model, yielding T G = 8.7^j j Gyr. 

Our two age estimates were combined using the maximum 
likelihood method, resulting in FINAL T c = (8.3 + 1.8) Gyr. 
This result is the first Galactic disk age determined via Th/Eu 
nucleocosmochronology, and is compatible with the most re- 
cent white dwarf ages determined via cooling sequence cal- 
culations, which indicate a low age (< 10 Gyr) for the disk 
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Determination of the Galactic disk age via [Th/Eu] nucle- 
ocosmochronology was found to be very sensitive to observa- 
tional uncertainties. The GDCE model is very insensitive to the 
choice of disk age, and the curves derived from the model are 
very close when compared to the abundance uncertainties. This 
leads to an age uncertainty as high as 1 .9 Gyr relative to the 
abundance ratio uncertainties alone, not taking into account the 
uncertainties intrinsic to the GDCE model itself, which are very 
difficult to evaluate. The analysis performed with production 
ratios and solar abundances taken from the literature presents 
an even higher uncertainty. However, considering that we man- 
aged to reduce the Eu and Th abundance scatters significantly, 
when compared to the best data currently available in the lit- 
erature, we believe that additional improvements may raise the 
precision of the analysis even further. Among the possible en- 
hancements, we can cite the observation of higher resolution, 
higher S/N ratio spectra (preferentially with large sized tele- 
scopes), determination of higher precision atomic parameters 
(e.g., central wavelengths of absorption lines), and identifica- 
tion of yet unknown line blends, like those that forced us to 
include artificial Fe lines in the spectral syntheses. If part of 
the data scatter is real, and not observational, it will not be pos- 
sible to reduce it indefinitely. This may be true if the Th/Eu 
production ratio is not constant throughout the Galactic evolu- 
tion, as abundance analyses of r-process elements in CS 3 1 082- 
001 seem to indicate. Future advancements in r-process nucle- 
osynthesis models may help solve the issues of universality and 
constancy of the production ratio, especially when the sites of 
production by this process are finally identified. 
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